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Q ' We show that the Kobayashi-Maskawa (KM) parametrization of the 3x3 lepton flavor mixing 

O ' matrix is a useful language to describe the phenomenology of neutrino oscillations. In particular, 

, it provides us with a convenient way to link the genuine flavor mixing parameters {9i, 62, ^3 and 

^km) ^° their effective counterparts in matter {61, O^, 0^ and (5km)- rediscover the Toshev- 
like relation sin Jkm sin 2^2 = sin (J^m the KM parametrization. We make reasonable 

I analytical approximations to the exact relations between the genuine and matter-corrected flavor 

mixing parameters in two different experimental scenarios: (a) the neutrino beam energy E is above 
C'(l) GeV and (b) E is below 0(1) GeV. As an example, the probability of i'^ oscillations 
^ ■ and CP-violating effects are calculated for the upcoming NO^A and Hyper-K experiments. 
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1 Introduction 



Just like quark flavor mixing, lepton flavor mixing has been observed in a number of neutrino oscillation 
experiments [1]. The phenomenon of lepton flavor mixing is described by a 3 x 3 unitary matrix V, 
the so-called Maki-Nakawaga-Sakata-Pontecorvo (MNSP) matrix [2]. There are totally 9 different 
parametrizations of V in terms of the rotation angles and phase angles [3]. Among them, the one 
advocated by the Particle Data Group (PDG) [1] is most popular in accounting for current neutrino 
oscillation data, and the Fritzsch-Xing (FX) parametrization has a particular merit in describing the 
running behaviors of neutrino masses and flavor mixing parameters from one energy scale to another 
by means of the one- loop renormalization-group equations [1]. Is the original Kobayashi-Maskawa 
(KM) parametrization |5] advantageous to the description of neutrino phenomenology? We shall give 
an affirmative answer to this question in the present work. 

In fact, it has recently been noticed that the KM parametrization of the 3x3 quark flavor mixing 
matrix is very useful to link current experimental data to the unitarity triangles — the so-called 
"Unitarity boomerang" [6j: ~ a^rj, ~ 90°, where is the CP-violating phase in the KM 

parametrization and a^T inner angles of the KM unitarity triangles. A similar relationship 

was found earlier in the FX parametrization [7] . As the structure of the KM parametrization is partly 
analogous to that of the FX parametrization and partly analogous to that of the PDG parametrization, 
we naturally expect that it should also be useful to describe the salient features of neutrino oscillations 
in vacuum and in matter. In other words, we expect that the KM parametrization can provide us 
with a simple and convenient link between its parameters and the observable quantities of neutrino 
oscillations. The main purpose of this paper is just to demonstrate our expectation and offer an 
alternative description of lepton flavor mixing for both the phenomenology of neutrino oscillations 
and the building of neutrino mass models. Needless to say, a convenient parametrization is sometimes 
possible to make the underlying physics more transparent. 

The remaining parts of this work are organized as follows: 

(1) In section 2 we shall first establish the explicit relations between the flavor mixing parameters in 
the KM parametrization {O^, 0^ and ^km) those in the PDG parametrization {9^2? ^13! ^23 

5). Because the Majorana CP-violating phases have nothing to do with neutrino oscillations, they will 
not be taken into account in this work. We find that 62 = ^23 '^km ~ ^ hold exactly if the condition 
cos 6 = sin cot 9^2 cot 2^23 (o^' equivalently cos S^^j^ = tan 63 cos 6-^ cot 262) is satisfied. Because both 
9^^ and 6*3 are expected to be small, the above condition seems to hint at (5 ~ ~ ±90°. 

(2) Section 3 is devoted to a detailed calculation of the relations between the genuine KM flavor 
mixing parameters in vacuum {9i, 92, 9^ and S-^j^) and their effective counterparts in matter (9i, 
^2) ^3 and (^km)' given neutrino oscillations in a constant terrestrial matter profile. The so-called 
Toshev relation sin 5 sin 2^33 = sin 5 sin 2^23 [S] in the PDG parametrization is rediscovered in the 
KM parametrization: sin (5km sin 2^2 = sin (5j^j^ sin 2^2 • This interesting result means that the KM 
parametrization is definitely useful and convenient to describe the phenomenology of neutrino oscilla- 
tions. To be more explicit, we make analytical approximations for the relations between the genuine 
and matter-corrected parameters of lepton flavor mixing in two different experimental scenarios: (a) 
the neutrino beam energy E is above 0(1) GeV; and (b) E is below 0(1) GeV. The accuracy of each 
approximation is examined by comparing its results with exact numerical calculations. Such analytical 
results are phenomenologically useful, just like those obtained previously in the PDG parametrization. 

(3) For illustration, we consider u — )■ oscillations in section 4 and calculate the oscillation 
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probability by means of the KM parametrization. Our expressions are simple and instructive, and they 
can be used to analyze the upcoming data from T2K [9], NOi/A [TO] and Hyper-K [11] experiments. 
(4) A brief summary of this work, together with some concluding remarks, is given in section 5. 



2 Comparison between the KM and PDG parametrizations 

In the framework of 3-generation leptons, the KM and PDG parametrizations of the MNSP matrix 
are given by 



V, 
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respectively, where Sj = sin^j, q = cos^j, = sin^j^ and = cos^j^ (for i = 1,2,3 and ij = 
12,23,13). Comparing the KM parametrization with the PDG parametrization, one can derive the 
expressions of and terms of 0j • and b: 
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With the help of Eq. dH 
properties of 0j and ^km- 



1(1 — sin 0;^3 cot cot ^23^"^) (1 + ^in ^13 cot 9^2, tan 023^"^) I 
and the current experimental oscillation data [12j, we find the following 



^1 ~ ^12 is ^ good approximation because of ^^3 < 12° as constrained by the present experimental 
data 



^2 depends on not only ^23 but also Q^z ^'^d b. 62 ~ ^23 ^ good approximation if 5 is near 
±90°, but not good if 5 is near 0° or ±180°. If one requires 62 = 6*23 to hold exactly and O^^, / 0) 
a constraint equation in the PDG parametrization must be satisfied: 



cos 5 = sin cot 612 cot 2^23 



(3) 



When calculating the effective flavor mixing parameters and neutrino oscillation probabilities 
in matter, one of ten takes ^^3 ^ expansion parameter to make analytical approximations 
because of its small value {e.g.^ P^)- Since tan03 = tan0^3 csc0]^2 holds, we can also treat 9^ as 
a small parameter to do series expansions. 

The experimental data yield ^23 ~ 45° [12], leading to sin Jj^j^ w sin 5 ± 0(513) from Eq. 
So b-^yi ~ (5 is a good approximation. Under the condition of Eq. ([3]), b-^yi = <5 holds exactly. 
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We remark that Eq. ([3]) is a useful constraint that leads exactly to both 02 = ^23 ^km ~ ^- 
be re-expressed in the KM parametrization: 

cos (5km = ^3 cos 9i cot 202 . (4) 

Moreover, one can see from Eqs. ^ and @ that 62 = ^23 is equivalent to 5-^^ = 5. One special case 
is the /i-r flavor symmetry with 62 = ^23 = ^5° and maximal CP violation with 5^^^ = 6 = ±90° |14j . 
We shall consider this possibility in the next section when discussing matter effects. 

We obtain the best-fit values and the la, 2a and 3a ranges of mixing angles 9i, 62 and ^3 for both 
the normal hierarchy (NH) and inverted hierarchy (IH) of neutrino masses in Table 1, according to 
a global analysis of current neutrino data presented in Ref. [12]. Some details of this table should 
be mentioned. First, we do not list any possible values of the CP phase 6^^^ because 6^^^ is not well 
constrained by current experiments. Ref. [12] gives very loose bounds of 6, with 5 = —110° (—74°) 
for the NH (IH) at the best fit and -227° to +7° (-200° to +43°) for the NH (IH) in the la range, 
without constraints in the 2a or 3a ranges. In this case, one may assume the bounds of 6 as the 
bounds of S^j^ for ^j^j^ ~ 6. Second, because we have little knowledge about 6, 62 cannot be well 
restricted from Eq. ([2]). Thus we assume 6 at its best-fit value when we calculate the range of O2 in 
Table 1. For a full understanding of the possible values of 62 and their dependence on 6, one may see 
Fig. 1, where we take the NH as an example to show how 62 changes with 5. 



3 The KM flavor mixing parameters in matter 
3.1 General formalism 

To understand the phenomenon of neutrino oscillations in the long-baseline experiments, it is neces- 
sary to analyze neutrino mixing in matter. In this section, we calculate the effective flavor mixing 
parameters in a constant terrestrial matter profile {0^, 62, 6^ and and study their relations with 
the genuine flavor mixing parameters in vacuum {6^^, 62, 9-^ and S^^^). 

In the flavor basis \i^{t)) = (|i^e(^)) A'^^ii'^)) )l'^r(^)))^> evolution of neutrinos in matter is 
described by a Schrodinger-like equation: 



with the effective Hamiltonian 
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(6) 



Here V is the MNSP matrix, a = Am2i/ Am^i is the mass hierarchy parameter with the mass-squared 

J, and A = 2Ea/ Am'^i is a dimensionless variable arising from the matter- 



differences Am?- 

'J 



induced effective potential a = \/2GpN^ [15j. is the number density of electrons in matter, and it 
can be taken to be half of the number density of nucleons of the Earth. In this case, A is given by 



A ^ 0.085 



2.5 X 10-3 eV" 
Aml-^ 



1 GeVy V2.8g/cm3 



(7) 
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where is the mass density along the path of neutrinos. In most long-basehne neutrino osciUation 
experiments, is approximately a constant [9l[10l[TT]. Eq. ^ holds for neutrinos. When considering 
the evolution of antineutrinos, we have to perform the replacements V ^ V* and a —a in the 
effective Hamiltonian. 

In view of Eq. ([6|) , one may define the effective MNSP matrix V through 



H 



2E 




(8) 



where Aj are the eigenvalues of the matrix in the square bracket of Eq. ([6]). The effective matter- 
corrected mass-squared differences are written as Am^i = Am|]^ (A2 — A^) and Am|^ = Am|^ (A3 — A^). 

To describe the MNSP matrices V and V in the KM parametrization, we need three real rotation 
matrices and one diagonal phase matrix O^, O2, O3 and 




in which Oi is a rotation matrix in the (1,2) plane and O2 ^"^^ rotation matrices in the (2,3) 

plane. Thus V and V can be parametrized as 
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V = 0,UsO,0, , 

V = d^Usd.d, , 

with the flavor mixing parameters in vacuum (6^, 62, 9^ and S^-^ 
respectively. 

In the representation of the rotation matrices and phase matrix, Eq 
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and in matter {0-^, 02, 9^ and (^km)' 
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Being a real symmetric matrix, M can be diagonalized by an orthogonal matrix V = 020iO^ with 
three rotation angles 6^, 62 and 6^. Thus, H is diagonalized by V' = 02Ug020iO^. After a phase 
transformation, one can derive V = U'V' = U' 0211^020^0^ with an additional unphysical phase 
matrix U'. Comparing it with Eq. ([9|), we find the following relations: 



^3 ~ ^3 ' ^2 — '^2^2 + *2^2 + 202*2^2^2 cos Jj^M 



2a2 



'1 > 



and 



sin 6km sin 2^2 = sin d-^u s™ 2^2 > 



(12) 



(13) 



where §2 = sin ^2 , and C2 = cos ^2 • 
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Eq. (jl3p gives a simple relationship between the genuine and matter-corrected flavor mixing 
parameters, which is similar to the Toshev relation [8] in the PDG parametrization 

sin J sin 2^23 = sin (5 sin 2^23 > (14) 

in which the tildes always stand for the parameters in matter. We point out that such a similar 
relationship exists in not only the PDG parametrization, but also the KM parametrization and another 
parametrization denoted as P7 in Ref. [3j Q. Generally, we call all of them the Toshev relations. 

We end this part with a direct application of this relation. Given the /i-r symmetry with 62 = 45° 
and maximal CP violation with S-^y^ = ±90° in vacuum, both sides of the Toshev relation equal 
±1, then |sin(5KM| = |sin202l = 1, leading to O2 = 45° and 5km = ±90°. Thus, we have proved 
that the /i-r symmetry and maximal CP violation keep unchanged when matter effects are taken into 
account. Together with the PDG and P7 parametrizations, we rewrite the overall invariance of the 
jjL-T symmetry and maximal CP violation in three parametrizations as follows: 

^23 = 45°, (5 = ±90° ^ ^23 = 45°, (5 = ±90° 

^62 = 45°, 5km = ±90° ^ '62 = 45°, 5km = ±90° 

^ e'2^ = 45°, 5' = ±90° ^ e'2Z = 45°, 5' = ±90° . (15) 

Ref. |16j provides another proof to this claim in the PDG parametrization, where a special basis of 
the neutrino fields is taken. 



3.2 Analytical approximations 

Our approximate formulation is based on two premises: 

• From Table 1, one can derive s| ^ 0.04 (0.05) for the NH (IH). So s| ~ \a\ ~ 0.03 provides a 
reliable basis for our analytical approximation. 

• Given small a, two scenarios should be considered separately: (a) E is above 0(1) GeV; and 
(b) E is below 0(1) GeV. 

The reason for distinguishing between scenarios (a) and (b) is simple. In scenario (b), where \a\ > \A\ 
holds, we have to regard ^ as a small parameter; but in scenario (a) with \a\ <^ \A\, we do not have 
to do so. 

^The P7 parametrization of the MNSP matrix is given by 



V 



C12C13 S12 -C12S13 

"■5'i2c'i3C23 + s'i3S23e c'i2C23 •5'l2s'l3C23 + c'i3S23e 



\ s'l2c'l3S23 + s'i3C23e — c'i2S23 ~'S'l2s'i3S23 + c'lsC^^ 



where s'ij = smO'ij and c'ij = cos O'ij. In this parametrization, the Toshev relation is expressed as sin 5' sin 2^23 ~ 
sin 5' sin 26'23 . The reason for the Toshev relation being satisfied is that the first rotation matrix on the right-hand side 
of the MNSP matrix is a (2,3) rotation matrix in all of these three parametrizations, and it commutes with the effective 
potential. 
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Scenario (a) 

Our first step is to diagonalize M in Eq. (Illj) . The symmetric matrix M can be decomposed into two 
terms in series of a: 



(16) 



The first term M(°) is a singular matrix and can be strictly diagonalized. The second term aM^^^ is 
a perturbation and can be diagonalized in series of a order by order. In general, the eigenvalues and 
eigenvectors of M are expressed as 



(0) , (1) , 



(17) 



respectively. And the orthogonal matrix V is written as V = {vi,V2,v^). Here one must pay attention 
to the detail of how to derive a proper order of Aj (or v^). In order to derive the proper order, we have 
to consider the cases of A < 1 and A > 1 separately. For ^4 < 1, we formulate the series expansion in 
a and derive the eigenvalues and eigenvectors to the first order. The eigenvalues are given by 
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The eigenvectors are too complicated to be listed here. 
Up to 0(si), Eq. (HID can be simplified to 
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The orthogonal matrix V turns out to be 
/ 
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whose three rotation angles are given by cot^^^ = ac^s^/A, tan ^2 = ~'^i^3 tan^g = SiS^/{l — A). 
Substituting the angles into Eqs. (fT2]) and (fT3]) . we obtain the approximate expressions of the matter- 
corrected flavor mixing parameters: 

a sin 26-, 



sint/2 = sint/21 



1 — 2e cos Sj^yi cot 62 + cot^ 6^ 



1 + 



tan^q — — 

I- A 



(l + e2) sin ,5km 



sin 5km = ' ^^"^ = , (22) 

1 + 2e cos (Jj^j^ tan 6*2 + tan^ 6*2) (l — 2e cos 6-^^ cot ^2 + cot^ ^2) 

in which e = c^^Sg. One can see that cot^^^ and tan^g are suppressed by a or S3, leading to 9i ~ 90° 
and 6*3 ~ 0. The most interesting result comes from 62 and ^km- Fhst, 62 ~ O2 and are 
good approximations for e is suppressed by S3. In other words, the matter effects on 02 and 6^^^ are 
small. Second, such small matter effects are independent of A in the approximation made above. This 
is a peculiar feature of the KM parametrization of V. 

Note that Eq. ()22p only holds for neutrinos with A < 1. For neutrinos with A > 1, the replacement 
of tan ^3 ^ cot ^3 in Eq. (f22l) should be made. 

For antineutrinos with —A<1, the approximate expressions are given by 

~ sin 01 sin^-n 



l + A ' 



sm tir, = cos ( 



'1 + 2e cos S-^yi tan 62 + tan^ 62 



l + e2 



~ a(l + A)cos9, 

cot Un = . . _ 

A sm 6*3 



(1 + e^) sin^T^.^ ^ ^ 

sin 5km = ^ ^ ™ — , (23) 

(1 + 2e cos (5j^]y[ tan 82 + tan^ 62) (l — 2e cos (5j^f^ cot 62 + cot^ Q2) 

in which e = C]^S3[1 — a(l + ^)/(As|)]. Eq. (p3|) is different from Eq. (j22]) because the order of 
and A2 has been exchanged. Given small ^3, ~ 0, ^3 ~ 90° and 62 ~ 90° ~ ^2 roughly hold. As e is 
^-dependent, 6*2 and obviously rely on A. 

For antineutrinos with —A > 1, the replacement of tan^^ =^ cot^^ in Eq. (|23p should be made. 

The magnitude of the intrinsic CP violation in neutrino oscillations depends only upon the Jarlskog 
invariant [17] . It is calculated via J' = Im (^1^2^*2^*1)- parametrization, one has J' = 

(1/8) sin 9^ sin 29^ sin 2^2 sin 2^3 sin^j^j^ in vacuum and J7 = (1/8) sin 9i sin 29^ sin 2^2 sin 2^3 sin(5j^j^ in 
matter. They satisfy the Naumov identity J^A?7i2iAm|]^Am32 = J'Am|]^Am|]^Am|2 [I8l[l9]. One 
can derive J' from Eq. (fT8]l as 

^ " AC(l-sin2 0iSin2 03)^ ^^^^ 

for neutrinos. Suppressed by a, the Jarlskog invariant in matter is much smaller than that in vacuum. 
J obtains its relative maximum J'^ ^ aj / {2s^s^) = (l/8)asin20;^ sin202 cos03sin(5j^]yj at A 1, 
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where C is near its minimal value. The series expansion in Sg gives the leading-order result of J 

a 



J 



:J ■ 



(25) 



A{l-Ar 

It can also be obtained from Eq. ()22p . Eq. ()25p does not hold for j4 ~ 1. To derive the correct result 
of J for antineutrinos, one has to make the replacements A —A and J — J . 

Scenario (b) 

Considering about \A\ < jaj, we rewrite Eq. ([TB]) as 

M = M'^°^ + AM'^^'^ + aM^^) , (26) 

where M'^^^ + AM''^^^ = M^^\ M'^^"^ has two degenerate eigenvalues and one eigenvalue 1. To the 
first order of a, the eigenvalues are expressed as 

\^ = ]^[A{l-slsl) + a-D] , 
1 



\2 = l^[A{l-slsl) + a + D] , 
A3 = 1 + VIS1C3S3 , 



(27) 



in which D 



[a (1 — slsf) + a] — AaAcf. M is diagonalized by 
/ cos -i? — sin 'd 



V = 0,0, 
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with one possibly large rotation angle •& at the zeroth order in A and two small perturbation coefficients 
/i^3 and /123 at the first order. A strict calculation yields 
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23 
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sint? 



1 A cos 26*1 - a 
- H 

2 2D 



= AsiS^{ci sin I? + SiC^ cos '&) , 
/123 = AsiS3{ci cos ?? — S1C3 sin -i?) . 



(29) 



Finally, by using Eqs. (jl2p and ()13p and taking account of ,83 ~ a, we obtain the flavor mixing 
parameters in matter: 
a 

sin 29i = — sin 26i , 



sin( 



tan 0o 



sin 6km 



sm I 



1 -2k cos cot ^2 + cot^ I 



1 + 



smfi smf 



(1 + K^) sin (5; 



KM 



1 + 2k cos (5j^]y[ tan 0^, + tan^ 0^,) (l — 2k cos (5^^/, cot 0^ + k? cot^ 



(30) 



KM ' 



where k = S3 sin esc (0^^ + t?). If one sets ^ — t- 0, then D — t- a, — t- and all the parameters in 
matter return to those in vacuum. For a nonzero A, 0^ may have a remarkable deviation from 0^. For 
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example, if A = q cos 20^, then D ~ Q!sin2^-|^ and sin 20^ ~ 1, leading to 9^ ~ 45°. But if A = a, then 
D ~ 2as-^ and sinlO-^ w cos0^, leading to 9^ w 45° + 0^/2. 9^ remains small though it receives some 
corrections. 6*2 and have small vl-dependent corrections suppressed by Sg. 

Eq. (j30p holds for neutrinos with the NH. One can replace i9 with 90° — i?, — and "(9 — 90° but 
keep Eq. (j30p unchanged to derive the correct results for neutrinos with the IH, antineutrinos with 
the NH and antineutrinos with the IH, respectively. 

The Jarlskog invariant reads 

'^ = ^'^ (31) 

for the NH. Since ~ A ~ a, ^7 is of the same order as J. Therefore, the CP violation for E < 0(1) 
GeV is not suppressed significantly by matter effects. J gets its relative maximum ~ esc 29^ — 
(1/8) sin 9^ sin 26*2 sin 26*3 sin (Jj^j^ at j4 ~ a cos 29 -^^j which is even greater than Jj. One can obtain the 
results in the IH case by changing a to —a. 

In Table 2 we show how to obtain the expressions of the matter-induced flavor mixing parameters 
in all the possible cases we have discussed in both scenarios. 

3.3 Numerical Analysis 

In Fig. 2 we make a comparison of the analytical and numerical results for the matter-corrected flavor 
mixing parameters. The analytical expressions in different intervals of A are given in Table 2. We 
have utilized Eq. with = 2.8 g/cm^ and replaced A with the neutrino energy E. The values of 
the other input parameters are taken to be Amgi = 7.59 x 10~^ eV^, Amg^ = 2.5 (—2.4) x 10~^ eV^, 
9^ = 34.5° (34.6°), 9^ = 49.3° (43.2°), ^3 = 11.6° (12.9°) and J^m = "HO" (-^4°) for the NH (IH), 
according to the best-fit values in Ref. [12] and Table 1. Some comments and discussions are in order. 

• In most cases, the analytical results in scenario (a) fit well with their numerical results for 
E > 0(1) GeV and the analytical results in scenario (b) fit well with their numerical results 
for E < 0(1) GeV. So |q;| <C |v1| corresponds to E > 0(1) GeV, and |q| > \A\ corresponds to 
E < 0(1) GeV. 

• The expression of 9^ for neutrinos with the NH and that of 9^ for antineutrinos with the IH are 
not good approximations in the interval 9 GeV < E < 15 GeV. Both 9^ and 9^ increase rapidly 
from a small angle to near 90° when E is running from 9 GeV to 15 GeV. 9^ reaches 45° at 

~ 12 GeV (or yl ~ 1) for neutrinos with the NH, while 9^ reaches 45° at the same energy for 
antineutrinos with the IH. 

• The analytical approximations of ^2 ^km good agreement with their numerical results. 
For neutrinos with E > 0(1) GeV, the numerical results confirm that ^2 ^''^^ <^km h^-ve small 
deviations from the vacuum parameters ^2 ^^'^ <^km ^^'^ nearly independent of the neutrino 
beam energy E. 

In Fig. 3 we compare the analytical results of the Jarlskog invariant J in Eqs. ()24p . (j25p and 
(131 p with the numerical results. The input values are taken the same as in Fig. 2. The analytical 
approximations fit well with the numerical results, except for Eq. ()25p in the interval 9 GeV < E < 
15 GeV. The numerical results show that the relative maximum of the Jarlskog invariant ~ —0.003 
takes place at ~ 12 GeV (or ^ ~ 1) for neutrinos with the NH and antineutrinos with the IH, and 



10 



~ —0.03 at ~ 0.1 GeV (or A ~ a cos 2^^) for neutrinos with both the NH and IH. They verify 
our analytical calculations of and in the above discussions. 

4 The probability of oscillations 

The probability of — t- oscillations in matter has been calculated in the PDG parametrization 
]yi\ [20t EH [22] but not yet in the KM parametrization. In this section, we calculate it in the KM 
parametrization. Both scenarios (a) and (b) will be considered. 
In vacuum, the oscillation probability of — )• is given by 

3 3 

P{y^^v:) = -4Y,Re{V^^Ve,V;^V:^s■mA,^-8JllsinA,^ , (32) 

i>j i>j 

in which A - = Amf-L/{'iE). In the KM parametrization, one has 

P {y^ —7- z^g) = (sin^ 9^ siv? 62 s\v? 26^ + J sin^ 6^ cos (5j^p^ + 2 sin^ 29-^^ cos^ 62 sin^ ^3) sin^ A32 
+J cos {Ar^i + (^km) ^32 ^^21 + sin^ 26*^ cos 62 cos^ 6*3 sin^ A21 
+2 sin^ 29i cos^ 62 sin^ ^3 cos A.^^ sin A32 sin , (33) 

where J = 8J / sin (5j^]y[ — sin sin 20]^ sin 2^2 sin 2^3. Replacing d^y^ with ~5kM' may obtain the 
expression of P (y^ — )• z^g) . Their difference is a measure of the intrinsic CP violation: 

AP {v^ ^ l/g) = P {v^ ^Ue)-P (P^ Ve) 

= — 2 Jsin^KM sin A3;^ sin A32 sin A21 . (34) 

By replacing the parameters in vacuum with those in matter, one can achieve a similar expression 
of the oscillation probability in matter P {v^ — )• from Eq. (]33p . In scenario (a), with the help of 
Eq. (I22p . P {v^ —7- fg) will be re-expressed in terms of the genuine mixing parameters and A. It is not 

difficult to derive the following approximate formula from Eq. ()33p by replacing A32 and Ag^ with 

sin(l-A)A3i , sin^Ao, . , 
— — — and a ^ , respectively: 

P (z^/. ^ ^e) = (sm^ ^1 sin^ 62 sin^ 2^3 + J sin^ ^3 cos 5km) ^'"^ [[Za)^^^ 

, , sinAAo, sin(l - ^)Ani 
+aJcos (A31 + 5km) ^ I- A 

siii"^ /4 /\ 

sin^ 201 cos^ 6^2 cos^ 6^3 ^^"^ ■ (35) 

Here A31 = A31, (1 — A)A.j^^ = A32 and ^A3i = A21 hold to the leading order in a and S3, and Eq. 



(|35p holds to the second order in a and the third order in S3. The second term in the first line is a 
term of S3. We keep it in view of the sizable value of ^3. The CP violation is included in the second 
line, which is suppressed by a. The term of in the third line is reserved for comparing with some 
former works in the PDG parametrization, in which the term of is often kept {e.g., [12 )■ But given 
the current experimental data, this term is meaningless unless we keep the terms of SI3 or asi3 |20j . 
The divergences for ^ — )• and 1 are absent because of the suppressions by sine functions sin AA and 
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sin(yl — 1)A. One can see from Eq. (j33p that the expression of the — t- oscillation probability in 
the KM parametrization is similar to that in the PDG parametrization. By replacing (5j^j^ with —li-^y^ 
and A with —A in Eq. (I35p . one can obtain the expression of P{h'^ — )• t^g) ^'^d then their difference 
AP {y^ — )• z^g) = — )• z^e) — P(^'^ — )• z^e)- '^-f' (^/i ~^ ^e) consists of both the intrinsic CP violation 
and matter- induced contribution and is sometimes dominated by the latter. 

In scenario (b), one can also make the similar replacements as in scenario (a) and obtain the 
approximate expression of P{v^ — s- in terms of the vacuum parameters and A. But because of the 
complexity of it is difficult for us to simplify the expression to a form as a simple function of A |21j . 
So we do not write out the expression here. We use another method proposed in Ref. [22] to obtain 
an alternative simple formula. It reads 

P {^t, ^^e) = P (j^M ^ ^e) + sin A31 (sin A31 - A31 cos A32) 

X (sin^ 9 1 s\T? 62 sin^ 2^3 + J sin^ ^3 cos 6km) ■ (36) 

With 

"^KM ^ ~'^KM ^ ^ '^^^ may arrive at P [v^ — )■ . Their difference is given by 

AP {v^ ^ u,) = AP [y^ ^ v^) + 4^ sin A31 (sin A31 - A31 cos A32) 

X (sin^ 0^ ■svT? 6*2 sin^ 26*3 + J sin^ 6*3 cos (5km) • (37) 

In Eq. ()37p . the intrinsic CP violation and the matter- induced contribution are separated to two 
different parts. The former is dominant if 6-^^ is not too small. 

Finally, we compare the analytical results of P {y^ — t- and AP {y^ — t- with the numerical 
results to show the validity of our approximations. We take the future long-baseline neutrino oscillation 
experiments NOi^A as an example in scenario (a) and Hyper-K as an example in scenario (b). The 
former will use the neutrino beams with E' ~ 2 GeV, which is compatible with scenario (a). And the 
latter will use the neutrino beams with E ~ 0.6 GeV, which is compatible with scenario (b). Different 
baseline lengths L and matter densities have been input in different experiments. L = 810 km 
and = 2.8 g/cm^ are taken in Fig. 4 for the NOz/A experiment [lOj, and L = 295 km and 

= 2.6 g/cm^ are taken in Fig. 5 for the Hyper-K experiment [11]. Other input parameters are 
taken the same as in Fig. 2. We see that the numerical results confirm the precision of the analytical 
approximations of P (f^ — t- z/g) and AP (z/^ — t- z/g) in both scenarios. 



5 Conclusions 

In this paper, we have used the KM parametrization to study the lepton flavor mixing and neutrino 
oscillations in matter. We re-discover the Toshev-like relation in the KM parametrization and prove 
that the p-T symmetry with 62 = 45° and maximal CP violation with (5j^j^ = ±90° keep unchanged 
when matter effects are taken into account. We have presented the approximate expressions of the 
matter-corrected flavor mixing parameters. Different methods have been chosen in two scenarios for 
the neutrino energy above 0{1) GeV and below 0(1) GeV. Finally, we have calculated the probability 
of — >• fg oscillations as an application in both scenarios. Below we compare the main features of 
the KM parametrization with those of the PDG parametrization: 

• The genuine flavor mixing parameters in the KM parametrization 9^, 62 and S-^y^ are approx- 
imately equal to the corresponding parameters in the PDG parametrization 9 12-, ^23 ^'^'^ 
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respectively. Although ^3 is not close to ^^3, it is small enough to do the series expansion in 
both ^3 and Am^i/ Am"^^ by treating sin^ 6*3 ~ Am^i/ Am^i. 

• For neutrinos with E above 0(1) GeV, the corrections to the rotation angle 62 and the CP phase 

induced by matter effects are small and nearly independent of the matter density and the 
neutrino energy. This is a salient feature of the KM parametrization. 

• The analytical expressions of the oscillation probabilities in Eqs. (I33p and (I35p are similar to 
the corresponding expressions in the PDG parametrization. 

It is well known that a good parametrization brings much convenience to the description of physical 
quantities. In this paper, we have explored the KM parametrization to study the neutrino phenomenol- 
ogy- 
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Parameter 




02 {S at its best fit) 


^3 


Best fit 


34.5° 


49.3° 


11.6° 




(34.6°) 


(43.2°) 


(12.9°) 


la range 


33.4° - 35.8° 


44.8° - 53.0° 


9.4° - 14.0° 




(33.5° - 36.0°) 


(39.8° - 45.8°) 


(10.4° - 15.3°) 


2(7 range 


32.1° - 37.4° 


41.9° - 55.0° 


6.8° - 16.0° 




(32.2° - 37.5°) 


(38.4° - 47.1°) 


(7.6° - 16.8°) 


3(7 range 


31.4° - 38.2° 


39.7° - 56.8° 


3.5° - 17.6° 




(31.4° - 38.4°) 


(37.8° - 48.2°) 


(3.5° - 18.6°) 



Table 1: The mixing angles in the KM parametrization translated from the results obtained in the 
PDG parametrization [12]. The upper (lower) row corresponds to the normal hierarchy (inverted 
hierarchy) of neutrino masses. The best-fit value 5 = —110° (—74°) have been taken for simplicity. 
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Table 2: The replacements for obtaining the proper expressions of the effective matter-corrected flavor 
mixing parameters in different cases. 

90° I ' 1 ' 1 ' 1 ' 1 ' 1 ' 
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6 

Figure 1: The dependence of 62 on 6 according to Eq. where the best-fit values and the la, 2a 
and 3a ranges of 612, ^13 and have been input. 
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Figure 2: The analytical (dashed line) and numerical (solid line) results of the matter-corrected flavor 
mixing parameters. The genuine flavor mixing parameters are taken as the best-fit values in Ref. 
[I2] and Table 1: Amj^ = 7.59 x lO'^ eV^, Am|i = 2.5 (-2.4) x IQ-^ eV^, 6^ = 34.5° (34.6°), 
6*3 = lf.6° (12.9°), 02 = 49.3°(43.2°) and = -110° (-74°) for the NH (IH). The mass density of 
matter is assumed to be = 2.8 g/cm^. 
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Figure 3: The analytical results of Eqs. ([25]) and (f3T]) (dashed line), Eq. (p^ (circle curve) and 
numerical results (solid line) of the Jarlskog invariant. All the input parameters are taken the same 
as in Fig. 2. 
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Figure 4: The analytical (dashed line) and numerical (solid line) results of the oscillation probability 
Pi^^fi ~^ ^e) and the difference AP(f^ — )• v^) with a length of baseline L = 810 km and matter density 
Pjjj = 2.8 g/cm^ in the NOi^A experiment |10) . The genuine flavor mixing parameters (0^, 62, 0^ and 
'^km) ^^"^ mass-squired differences (Am^]^ and Am^]^) are taken the same as in Fig. 2. 
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Figure 5: The analytical (dashed line) and numerical (solid line) results of the oscillation probability 
P{v^ — )• Ve) and the difference AP(z/^ — )• Uf.) in the Hyper-K experiment [11]. Compared with Fig. 4, 
a length of baseline L = 295 km and matter density p^^ = 2.6 g/cm^ are taken. The genuine flavor 
mixing parameters and mass-squired differences are taken the same as in Fig. 2. 
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